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Abstract — The  HF  ocean  radar  was  developed  in  Japan  from 
1987  by  Okinawa  Radio  Observatory,  Communications  Research 
Laboratory.  We  measured  ocean  surface  currents  and  developed 
the  method  to  estimate  wave  spectra.  The  observations  of 
ocean  surface  currents  were  conducted  in  various  areas.  The 
observation  in  the  east  of  Okinawa,  where  mesoscale  eddies  are 
dominant,  is  presented.  We  can  observe  the  convergent  zone 
associated  with  the  front  of  a  mesoscale  eddy.  Furthermore, 
the  short-wave  directions  and  spread  parameters  were  estimated 
from  first-order  scattering.  We  also  developed  the  method  to 
estimate  ocean  wave  spectra  from  Doppler  spectra.  We  extended 
the  theory  of  HF  radio  wave  scattering  from  the  sea  surface.  It 
is  possible  to  explain  the  underestimation  of  waveheights,  which 
was  previously  pointed  out. 

I.  Introduction 

HF  ocean  radar  is  used  to  measure  ocean  surface  currents 
and  waves  by  radiating  high-frequency  radio  waves  and  ana¬ 
lyzing  backscattered  signals  from  the  ocean.  This  technique 
is  based  on  the  backscattering  of  radio  waves  by  resonant 
surface  waves  of  one  half  the  incident  radar  wavelength.  This 
is  called  Bragg  scattering,  and  its  Doppler  frequency  (Bragg 
frequency)  is  determined  from  the  linear  dispersion  relation¬ 
ship  of  surface  waves.  However,  in  practice,  the  position  of 
first-order  scattering  of  the  measured  Doppler  spectrum  is 
somewhat  different  from  the  Bragg  frequency.  This  difference 
is  due  to  the  advection  of  the  surface  waves  by  ocean  currents. 
We  can,  therefore,  estimate  the  radial  component  of  the  ocean 
current  by  estimating  the  frequency  shift  of  the  Bragg  peaks. 
Furthermore,  the  two  first-order  Bragg  peaks  are  surrounded 
by  the  second-order  continuum,  which  is  called  second-order 
scattering.  All  of  the  wave  components  contribute  to  second- 
order  scattering,  and  we  can  estimate  ocean  wave  spectra  from 
second-order  scattering. 

The  development  of  HF  ocean  radar  in  Japan  began  from 
1987  in  Okinawa  Radio  Observatory  (Okinawa  Subtropical 
Environment  Remote-Sensing  Center),  Communications  Re¬ 
search  Faboratory.  Test  observations  were  conducted  only 
by  one  radar  from  1989.  We  developed  the  second  radar  in 
1991,  and  it  is  possible  to  measure  ocean  current  vectors.  We 
measured  ocean  surface  currents  and  developed  the  method  to 
estimate  wave  spectra.  The  objective  of  the  report  is  to  present 
past  and  present  studies  on  HF  radar  oceanography  in  Japan. 
The  topics  are  (1)  measurement  of  ocean  surface  currents  ([1]), 
(2)  short-wave  directional  properties  ([2]),  (3)  wave  directional 
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TABLE  I 

Characteristics  of  the  HF  ocean  radar 


spectrum  estimation  ([3]),  and  (4)  HF  radio  wave  scattering 
from  the  sea  ([4],  [5]). 

II.  Radar  system 

Table  1  is  the  specification  of  HF  radar  system.  The  radio 
frequency  is  24.5MHz  and  the  FMICW  (Frequency  Modulated 
Interrupted  Continuous  Wave)  radar  is  used.  The  phased  array 
system  is  used,  so  the  large  area  for  the  deployment  is 
necessary.  The  spatial  range  resolution  of  the  radar  is  1.5km. 
The  radar  beam  can  be  directed  from  — 45°to  •  15  in  7.5°step. 
The  beam  direction  is  controlled  by  the  hardware  in  real  time. 
It  takes  10  minutes  to  measured  one  direction,  therefore,  we 
repeated  the  measurements  every  2  hours  in  the  observation. 

III.  Observation  of  currents 
A.  Comparison  of  currents 

Figure  1  is  the  observation  area  by  HF  ocean  radar.  The 
radars  are  located  at  site  A  (26°7.19’  N,  127°45.78’  E)  and  B 
(26°  18.63'  N,  127°50.25’  E).  The  region  east  of  Okinawa  is 
affected  by  mesoscale  eddies.  The  HF  radar  observation  was 
conducted  from  15  April  1998  to  15  May  1998  in  the  east 
of  Okinawa  Island  (around  26°N,  128°E).  The  current  meter 
was  deployed  at  the  location  P  (Payao).  The  water  depth  of 
the  deployment  was  4  m.  The  wind  station  is  at  I  (Itokazu) 
in  Figure  1.  The  location  of  the  wave  measurement  is  at  Ky 
(Kyan)  in  Figure  1. 

We  compared  HF-ocean-radar-derived  currents  with  in-situ 
current  data  [  1  ] .  Table  II  shows  the  summary  of  the  com¬ 
parison.  The  number  of  samples  in  the  comparison  was  309. 
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convergent  zone  was  associate  with  a  temperature  front. 


Fig.  1.  Map  of  observation  area 


The  correlations  are  greater  than  0.85,  and  the  root-mean- 


Up  —  ur 

Vp  —  vr 

Correlation 

0.87 

0.89 

rms  difference  (cms-i) 

11.68 

10.47 

Slope 

1.015 

1.047 

Bias  (cms-1) 

-3.89 

0.043 

TABLE  II 

Comparison  between  HF  ocean  radar-derived  currents  (u,  : 

EAST- WEST  COMPONENT,  vr :  NORTH-SOUTH  COMPONENT)  AND  CURRENT 
METER-MEASURED  CURRENTS  AT  PAYAO  (Upt  EAST- WEST  COMPONENT, 
VP:  NORTH-SOUTH  COMPONENT). 


square  difference  between  the  sensors  is  about  11  cm  s-1.  The 
slopes  and  biases  of  the  regression  lines  to  the  scatter  plots  are 
respectively  about  1  and  —4  cm  s-1.  From  this  comparison  we 
conclude  that  radar-derived  currents  are  sufficiently  accurate 
for  our  purposes. 

The  HF  radar-observation  area  appears  to  be  affected  by 
the  mesoscale  eddy  front  in  the  observation  period  from 
altimetric  data  and  SST  (sea  surface  temperature)  data  by 
TRMM  (Tropical  Rainfall  Measuring  Mission)  microwave 
imager  (Figure  2b).  As  indicated  in  Figures  2,  prominent  con¬ 
vergent  zones  were  observed  off  Okinawa  Island.  For  example, 
the  southeastward  and  southwestward  flows  form  a  convergent 
zone  in  Figure  2a.  The  water  depth  of  the  convergent  zone  was 
greater  than  1000  m.  From  NOAA/AVHRR  imagery  data,  the 
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Fig.  2.  (a)  Daily  mean  currents,  (b)  daily  SST  data  by  TRMM  microwave 

imager. 


IV.  Short-wave  direction 


We  can  estimate  mean  directions  9m  and  spread  parameters 
s  at  the  Bragg  wavenumber.  They  are  estimated  from  the 
ratio  of  first-order  Doppler  spectra.  If  the  wave  directional 
distribution  is  written  as  D(6)  cx  cos2s (9/2),  where  0  is  the 
wave  direction  relative  to  the  mean  direction,  the  first-order 
ratio  is  written  as 


D(tt  -  9m) 
D(-9m) 


(1) 


We  can  estimate  9m  and  s  by  two  radars  from  (1).  The  short¬ 
wave  direction  is  an  indication  of  wind  direction.  We  compared 
short-wave  directions  with  wind  directions.  The  comparison 
is  shown  in  Figure  3.  The  agreement  is  good.  The  correlation 
was  about  0.96,  and  the  root-mean-square  difference  from  the 
linear  regression  line  was  about  31°.  The  difference  between 
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Scatter  diaqr  am 


Winds  at  Itokazu 


Fig.  3.  Scatter  diagram  between  wind  directions  with  respect  to  the  eastward 
direction  (counterclockwise  is  positive)  at  Itokazu  radar-estimated  short-wave 
directions  at  the  closest  grid  point  to  Itokazu. 


wind  directions  (I  in  Figure  1)  and  radar-estimated  short-wave 
directions  are  due  to  horizontal  variations  of  winds  and  a  time 
lag  of  short-wave  response  to  a  sudden  wind  shift.  If  the  wind 
held  is  steady,  the  rms  difference  is  small.  For  example,  a 
previous  study  showed  that  the  rms  difference  between  the 
radar  and  buoy  measurements  was  24° ,  when  high  quality  data 
were  used  [6], 

Figure  4  shows  radar-estimated  short-wave  directions  in 
25  April.  The  wind  direction  and  short-wave  directions  were 
northeastward  until  6  JST  (Figure  4(a)).  The  short-wave  direc¬ 
tions  in  the  northeastern  part  of  HF  ocean  radar  observation 
area  were  changed  to  southeastward  (Figure  4(b)).  We  can  see 
a  discontinuity  of  short-wave  directions  in  Figure  4(b).  This 
is  associate  with  a  passage  of  an  atmospheric  front. 

V.  Wave  spectrum  estimation 

We  developed  the  inversion  method  to  estimate  ocean  wave 
spectra  from  second-order  Doppler  spectra  [3],  The  integral 
equation  which  relates  Doppler  spectra  to  the  wave  spectrum 
is  discretized.  Thus  estimation  of  wave  spectra  is  converted  to 
the  optimal  problem.  We  demonstrated  from  the  twin  pair  sim¬ 
ulation  that  this  method  could  estimate  ocean  wave  directional 
spectra  by  two  radars  [3].  Figure  5  shows  an  example  of  wave 
spectrum  estimation  by  dual  radars.  It  is  possible  to  estimate 
frequency  spectra  even  if  the  Doppler  spectrum  is  available 
in  only  a  single  direction.  In  this  case,  however,  the  solution 
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Fig.  4.  Radar-estimated  short-wave  directions  in  25  April  1998  and  wind 
vector  at  Itokazu.  (a)  6  JST,  and  (b)  8  JST. 


of  the  directional  spectrum  tends  to  converge  to  a  spectrum 
that  is  symmetrical  to  the  beam  direction.  Even  if  the  wave 
spectrum  is  dominant  in  a  single  direction,  the  solution  may 
give  two  peaks  in  the  wavenumber  spectrum.  One  of  them 
is  the  true  peak  and  the  other  is  the  mirror  image  of  it  with 
respect  to  the  beam  direction.  It  is  necessary  to  use  two  beam 
directions  for  estimation  of  wave  directional  spectra. 

We  compared  buoy-measured  with  radar-estimated  ocean 
wave  spectra.  The  HF  radar  observation  was  conducted  in  the 
Japan  sea  from  5  to  20  March  1991,  and  the  radar  location 
was  at  Yura  (near  40°  N,  140°)  [7],  We  used  only  one  radar,  so 
we  can  only  estimate  frequency  spectra.  The  frequency  spectra 
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(a)  Retrieved  wave  directional  spectra 


Normalized  wave  frequency(coN) 
(b)  True  wave  directional  spectra 


Normalized  wave  frequency(«N) 

Fig.  5.  Comparison  of  wave  directional  spectra  as  in  the  dual-beam 
case  (crossing  angle  is  90°).  (a)  Retrieved  directional  spectra,  and  (b)  True 
directional  spectra  (thick  dotted  line=  2;  thick  dashed  line=  3;  thick  solid 
line=  4;  contour  interval^  0.5). 


were  integrated  and  were  converted  to  significant  waveheights 
by  calculating  the  square  root  of  the  integrated  values  and 
multiplying  by  the  factor  4.  Table  III  shows  the  summary  of 
the  comparison.  The  number  of  samples  in  the  comparison 
was  173.  The  slope  (a)  and  bias  (b)  of  the  regression  line 
Hr  =  a  1 1 1,  +  b,  where  //,.  and  Hb  are  respectively  radar- 
estimated  and  buoy-measured  significant  waveheights,  were 
estimated.  The  correlation  coefficients  are  more  than  0.8  for 
total  frequency  band,  and  are  large  enough  for  practical  use. 
However,  there  are  bias  errors  resulting  in  overestimation  of 
radar-estimated  waveheight.  The  bias  errors  are  primary  cre¬ 
ated  by  external  noise  in  Doppler  spectra.  The  noise  originates 
from  natural  (atmospheric)  and  man-made  (environmental) 
sources.  A  method  to  reduce  noise  needs  to  be  developed. 


frequency  band 

Correlation 

slope  (a) 

bias  (6,  m) 

0.03-0.10  Hz 

0.804 

1.058 

0.387 

0.10-0.125  Hz 

0.897 

1.023 

0.310 

0.125-0.155  Hz 

0.838 

0.858 

0.344 

0.16-0.20  Hz 

0.692 

0.808 

0.362 

0.20-0.25  Hz 

0.638 

0.754 

0.347 

Total 

0.858 

0.948 

0.718 

TABLE  III 

Comparison  of  wave  spectra. 


VI.  Radio  wave  scattering 

The  HF  radio  wave  scattering  theory,  which  is  the  basic 
equation  to  estimate  ocean  wave  spectra,  is  derived  from  the 
second  order  perturbation  method  of  electromagnetic  field. 
However,  we  can  expect  that  this  theory  does  not  apply  in  the 
case  of  high  sea  conditions.  The  radio  frequency  allocation 
is  quite  limited,  and  it  is  not  feasible  to  use  both  high-  and 
low-radio  frequency  according  to  the  sea  states.  That  is,  it  is 
not  feasible  for  us  to  change  radio  frequencies  according  to 
the  sea  states.  In  the  perturbation  theory,  the  scattered  electric 
field  is  written  as  E  =  E\  +  E2  +  E3  +  •  •  •,  and  the  radar 
cross  section  is  proportional  to  the  power  spectrum  of  the 
scattered  electric  field.  If  the  first-order  scattered  field  E\  is 
a  Gaussian  process  with  zero  mean,  the  power  is  written  as 
( \E\2)e  =  (|£i |2  +  \E2\ 2  +  +  E{E3  +  •  •  where  an 

asterisk  denotes  the  complex  conjugate  and  ■  -)e  denotes  the 
ensemble  average.  The  term  {\Ei\2)e  corresponds  to  the  first- 
order  radar  cross  section  and  the  term  {\E2\2)e  corresponds 
to  the  second-order  radar  cross  section.  The  terms  (E\E3)e 
(and  {El  E3)e)  are  of  the  same  order  as  the  second-order  radar 
cross  section.  These  terms  modify  the  first-order  radar  cross 
section,  and  they  affect  the  estimation  of  waveheight,  because 
the  second-order  radar  cross  section  is  normalized  by  first- 
order  radar  cross  section.  We  calculated  this  correction  term  of 
first-order  radar  cross  section  [4],  Figure  6  is  the  result  of  the 
ratio  of  the  correction  term  to  the  first-order  radar  cross  section 
as  a  function  of  dominant  wave  frequency  and  direction.  For 
example,  0. 1  Hz  of  the  dominant  wave  frequency  corresponds 
to  4  m  of  the  significant  waveheight,  the  ratio  is  about  0.2.  This 
means  that  if  the  correction  is  not  considered,  the  waveheight 
is  under-estimated  by  about  10  %.  The  under-estimation  of 
waveheight  is  consistent  with  [8]. 

It  is  for  higher-order  correction  of  the  first-order  Bragg 
echoes,  and  it  has  no  relation  with  the  shape  of  second- 
order  scatter  scattering.  If  we  see  Doppler  spectrum  at  high 
sea  states,  the  shape  of  the  Doppler  spectra  is  different  from 
predicted  Doppler  spectrum.  In  general,  observed  Doppler 
spectra  are  broaden  at  high  sea  states  [7], 

In  order  to  explain  this  discrepancy  between  measured  and 
predicted  Doppler  spectra,  we  calculated  second-order  Doppler 
spectra  by  considering  the  finite  illuminated  area  [5].  Figure  7 
is  the  result  of  the  comparison.  This  is  the  VHF  case  and 
the  radio  frequency  is  42  MHz.  In  these  cases,  the  prediction 
appears  to  be  improved.  But  there  are  still  discrepancies  as 
Figure  7a.  This  prediction  should  be  improved  further. 
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CONTOUR  INTERVAL  =  2.000E-01 


Fig.  6.  Third-order  correction  term  (corresponding  to  ( / ■.  i  IE)  /;  and 
(E*Es)  /  )  relative  to  the  first-order  radar  cross  section  as  a  function  of  peak 
wave  frequencies  and  peak  wave  directions 


VII.  Summary 

The  report  is  summarized  as  follows: 

(1)  The  agreement  between  measured  and  radar-estimated 
currents  is  good.  (2)  The  convergent  zone  associated  with 
mesoscale  eddy  is  detected  by  HF  ocean  radar.  (3)  The 
agreement  between  wind  directions  and  radar-estimated  short¬ 
wave  directions  is  good.  (4)  We  developed  the  inversion 
method  to  estimate  wave  spectra.  Comparisons  with  in-situ 
measurements  are  necessary  further.  (5)  We  extended  the  HF 
radio  wave  scattering  theory,  however,  the  further  improvement 
of  the  theory  to  explain  the  broadening  of  second-order  radar 
cross  section  is  necessary. 
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